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The ground-state correlations are investigated for an isotropic transverse XY chain which is 
constrained to carry either a current of magnetization (J^^) or a current of energy (J^). We find 
that the effect of J^' ^ on the large-distance decay of correlations is twofold: i) oscillations are 
introduced and ii) the amplitude of the power law decay increases with increasing current. The effect 
of energy current is more complex. Generically, correlations in current carrying states are found 
to decay faster than in the — states, contrary to expectations that correlations are increased 
by the presence of currents. However, increasing the current, one reaches a special line where the 
correlations become comparable to those of the = states. On this line, the symmetry of the 
ground state is enhanced and the transverse magnetization vanishes. Further increase of the current 
destroys the extra symmetry but the transverse magnetization remains at the high-symmetry, zero 
value. 
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I. INTRODUCTION 

A general feature of nonequiUbrium steady states is the 
presence of currents (fluxes) of some physical quantities 
such as energy, momentum, charge, etc. Thus the study 
of nonequiUbrium steady states is, in some sense, a study 
of the effects of currents imposed on the system either by 
boundary conditions and driving fields or by competing 
dynamical processes. An interesting and much investi- 
gated effect of currents is the rather dramatic change in 
correlations. Namely, short-range correlations appear to 
change into long-range, power-law ones as the currents 
are switched on [Q,^. This is not entirely surprising in 
case of a global current since some conserved quantity is 
carried fast compared to diffusion and, in the absence of 
detailed balance, this may generate long-range effective 
interactions and, as a consequence, long-range correla- 
tions may appear jsj. It is clear, however, that the general 
picture should be more complicated since large currents 
often lead to chaotic behavior which in turn results in 
weakened correlations Thus, we believe, the under- 
standing of the interplay of currents and correlations in 
nonequilibrium steady states is a rather interesting and 
important task. 

In order to achieve progress one tries to investigate 
simple models and, indeed, a large number of classical 
statistical models have been introduced for the study 
of nonequilibrium steady states Unfortunately, far 
from equilibrium, classical systems are not constrained 
by conditions like detailed balance and there is much ar- 
bitrariness in defining the dynamics. In order to avoid 
such arbitrariness, we have started to investigate quan- 
tum systems where the time evolution is defined with- 
out ambiguity by the usual rules of quantum mechanics. 

Nonequilibrium steady states in a quantum system 



may be investigated by imposing a current on the system 
and studying the properties of the ground state thus gen- 
erated. As an example, we studied the transverse Ising 
model 1^ and found that, in the presence of an energy 
current, the exponentially decaying two-spin correlations 
changed into power-law form thus supporting the notion 
that switching on currents increases correlations. In view 
of the lack of detailed knowledge of the interplay of cur- 
rents and correlations, here we probe the generality of the 
above result on the example of the isotropic transverse 
XY chain. In this system, we have a global conservation 
not only for the energy - as in case of the transverse Ising 
model - but also for the transverse magnetization and, 
consequently, the effect of both the energy- and magne- 
tization currents can be investigated. The XY chain is 
also interesting because it has power-law correlations al- 
ready in the equilibrium state (i.e. in the state without 
any current) and so one might expect the system to be 
more sensitive to the introduction of currents. Indeed, in 
case of the energy current we find a rather complex be- 
havior (including an increase of the ground-state symme- 
try at special values of the current) and our findings for 
the steady state correlations are at variance with those 
for the transverse Ising model. On the other hand, the 
changes we observe in the correlations due to a magneti- 
zation current allow for a straightforward interpretation 
in terms of increasing correlations due to presence of a 
current. 

The effects of magnetization current are treated ex- 
actly (Sec. II) while the correlations in the presence of an 
energy current are calculated by a combination of ana- 
lytical and exact numerical methods in Sec. III. Summary 
and final comments are contained in Sec. IV. 
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II. MAGNETIZATION CURRENT 

Our starting point is the Hamiltonian of the d 
isotropic XY model in a transverse field h: 
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where the spins are represented by Pauli spin matrices s" 
[a = X, z) at sites I = 1, 2, of a one-dimensional 
periodic chain (s^_|_]^ = s"). The transverse field, /i, 
is measured in units of the Ising coupling, J, which is 
set to J = 1 throughout this paper. This model can 
be solved exactly 0| since it can be transformed by a 
Jordan- Wigner transformation into a set of free fermions 
with wavenumber fc and of energy 



cos k — h. 



(2) 



In this model, not only the total energy but also the 
z component of the total magnetization = J2e 
conserved. As a result, one can write down a continuity 
equation for the local magnetization s|: 
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and this defines the magnetization current through the 
bonds: 



A macroscopic current can now be defined as 
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and one can recognize as the Hamiltonian of the 
Dzyaloshinskii-Moriya interaction jsj . Somewhat surpris- 
ingly, the same expression emerged as the energy current 
in the case of the transverse Ising model Q . 

Our aim is to find the lowest energy state among the 
states carrying a given current. Since [if"^^, J*^] = 0, 
the problem can be solved using the Lagrange multiplier 
method, i.e. we diagonalize the following Hamiltonian: 



(6) 



where A is a Lagrange multiplier. The ground state of 
can be considered as a current-carrying steady state 
of i?^^ at zero temperature. Note that, without loss of 
generality, we can assume h> and A > 0. 

The Hamiltonian H^^ is diagonalized using the same 
transformations which diagonalize H-^^ , and we get the 
following spectrum in the thermodynamic limit: 
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where (p = arctanA and an effective field h = hcosip 
has been introduced. One can see that the spectrum is 
similar to that of H^^ with the wavenumber shifted by 
ip. It should be mentioned here that the above result 
in not new. It is implicit in the Bcthe-ansatz solution 
of the anisotropic Heisenberg chain with twisted toroidal 
boundary conditions and appears in various forms in 
studies of the Dzyaloshinskii-Moriya interaction and 
of the associated Berry phase Q . 

It is remarkable that H^' can be transformed into 
(pi) using a unitary transformation: 
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This transformation rotates the £th spin around the z 
axis by angle iip and shifts the spectrum (|^) by the 
wavenumber (p. This shift is analogous to the phase shift 
experienced by electrons on a ring threaded by a constant 
magnetic flux ($ ^ Nip) 

We also note that the ferro- and the anti-ferromagnetic 
cases are equivalent in the sense that the canonical trans- 
formation (j^) with (p = n transforms them into each 
other. For a finite periodic chain, H^^ is transformed to 
, but with a twisted boundary condition {jp ^ it). 

As the transformation (H) does not change the z com- 
ponent of the spins, we find that 
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and the correlation function 
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have their equilibrium form but at a different field h — 
/i/vl + A^ (( ) denotes the expectation value in the 
ground state of H^). One can see that, for h > Vl + A^ 
{h > 1), the spins are parallel to the field and there is 
no current in the system. For h < vl + A^ {f^ < ^) the 
current has a simple form: 



A^/TTA2~^ 
41 + A2) 



(11) 



and one can observe that the maximum current, reached 
in the limit of A ^ cx3, is given by I/tt. 

The introduction of the transformation (^) simplifies 
the calculation of the correlations p^{r) — (sfsf^^) = 
/9^(r) in the ground state of H^^ : 



p'^ir) = cos(r^)(QV'o|s?s?+rlWo) 
+ sin{rip) (QVo I s?s|_^ jg7/>o> 



(12) 



where ipo is the ground state of while Qtpo is that 
of H^^ at a field h. Without any current, we have 
(s|S^_|_^) = and, furthermore, the r — > oo behavior 
of the p^{r) correlation function jl^ is given by 
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p^{r;f^^O)^{l-hY'^, (13) 

where C = e^l'^^'^l^ A'^ « 0.147 {A w 1.282 is the 
Glaisher's constant). Using (p^, one can obtain then 
the following simple form in the r — s- cx) limit 

p"(r) w (l - ^'^ ^ cosM. (14) 

Thus we find that the correlations decay by power 
law and they show oscillatory behavior in the current- 
carrying states. This is similar to what has been ob- 
served in the transverse Ising model but there are some 
differences. In the Ising case, an exponential decay of cor- 
relations changes into a power-law form as the current is 
switched on. In the XY case, on the other hand, one has 
power-law correlations already in the equilibrium state. 
The magnetization current does not change the power 
law and leaves its exponent intact as well. The increase 
of correlations appears as the increase in the amplitude 
of the power law (note that h ~ h cos tp < h). 



h < 1) and 7r/2 ± k\ (for A > 1), respectively. The 'crit- 
ical momenta' and kx are defined here such that they 
take values < kh, k\ < 7t/2 and one has kh — arcsin (h) 
and k\ — arccos(A~^). 

One can study the ground state as a function of h and 
A, but we are more interested in the physical quantities as 
functions of h and = (J^/N). Thus first we calculate 

f = < ^Vl-A-2 for kh >kxOTh,X>l 

i for A < 1, 

(18) 

and than express all the A dependences in terms of . 
We can then obtain a. h — phase diagram as shown 
on Fig.|l] where the phases - discussed below in more de- 
tail - are distinguished by symmetries of the regions of 
occupied states in the fc-space. 



III. ENERGY CURRENT 



A. Phase diagram 



Since the energy is a conserved quantity as well, one 
can investigate the effects of the energy current. The 
local energy (the contribution of the £th spin to the to- 
tal energy) satisfies a continuity equation with the local 
energy current jf, and its sum over £ (the total energy 
current) has the form: 

JV 

(15) 

One can easily show that [H^^ , J^] — 0, and diagonal- 
izing the Hamiltonian 
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one obtains the lowest energy eigenstates of H^^ in the 
presence of a given . 

Using the standard transformations to fcrmions again, 
the spectrum is obtained as: 



Afc = (— cos k — h){l — X sin k) , 



(17) 



and the modes with negative energy are occupied in the 
ground state of . Although the k ^ —k symmetry 
of the spectrum is broken for A 7^ 0, the ground state 
remains that of H^^ for A < 1 and, accordingly, no en- 
ergy current flows through the system. This rigidity of 
the ground state against A is a consequence of the fact 
that the fermionic spectrum of has a product form 
(p^), and the second factor is positive for A < 1. The first 
and second factors in A^ change sign at ±{tt/2 + k^) (for 



Let us begin the analysis of the phase diagram by first 
describing it in terms of the behavior of currents and of 
the transverse magnetization. As shown in Fig.|l|, there 
is a maximal current for every value of h 
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and no state exists above j 



E 

max ■ 




FIG. 1. Phase diagram of the ground state of the trans- 
verse XY model in presence of energy current. The black 
parts of the rectangles denote the wavenumbers of the oc- 
cupied fermionic modes (— vr < k < tt). The dashed line is a 
high-symmetry transition line between regions of (T) (M^ — 0) 
and @ (M" 7^ 0). There are no states above the maximal 
current line and region ® can be mapped onto the vertical 
dotted line {h = 1). 
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We can divide the h-j^ phase diagram into three re- 
gions. The only interesting areas are Q) and @ and their 
boundaries. In region @ , the ground state is the same 
along the ~ constant x h lines {h > 1) thus the h = 1 
line contains all the information about this region. Be- 
low, we restrict the discussion to the h < 1 part of the 
phase diagram with the understanding that the h = 1 
line represents region @ . 

As can be seen from (^5|), the energy current has two 
parts: the term containing h is proportional to the mag- 
netization current {—hJ^^) while the other term is the 
current of the interaction energy. The distinguishing fea- 
tures of regions @ and (g) are that the current of inter- 
action energy is zero while the transverse magneti- 
zation, , is nonzero in the ground state. 

For any fixed value of h, the magnetization decreases 
with increasing and becomes zero on the line 
kh = k\ corresponding to = h'^/ir. On this line, 
the magnetization-current part of saturates and, upon 
increasing , we enter region (T) where the interaction 
part of the current starts to flow. Another characteris- 
tic feature of region (T) is that — throughout this 
region. 

One tends to conclude at this point that the line 
separating regions (T) and @ is a line of second order, 
nonequilibrium phase transitions with being the or- 
der parameter. This notion is also supported by the facts 
that several quantities such as and p^{l) have a 

jump in their first derivatives when crossing the tran- 
sition line and, furthermore, that the correlations are 
enhanced (change from to r~^/^) on this line (see 



Sec. pi B| ). If this was a phase transition, however, it was 
certainly a transition not in the usual sense. The sym- 
metry of the ground state is the same on both sides of 
the transition line and the = result in region (T) 
is not a consequence of the up-down symmetry of the 
ground state. The magnetization is zero in (T) because 
the motion of the zeros of the dispersion relation ( |l^ ) 
conspire to keep the ground state at half-filling. We em- 
phasize, however, that the half-filling does not mean that 
the ground state has a symmetry with respect to global 
spin- flip sf —sf. 

It is interesting to note that the symmetry of the 
ground state is higher on the transition line than on ei- 
ther side of it. Indeed, on this 'high'-symmetry line, the 
ground state is symmetric with respect to rotation of the 
spins around the x axis by tt, followed by a spatial reflec- 
tion mapping site i to L -I- 1 — i. The Hamiltonian, , 
has no such symmetry and, off the transition line, the 
ground state doesn't have such symmetry either. Thus 
we can see here an example where the increase of current 
in a system leads, at a particular value of the current, to 
symmetry enhancement. The reason for increase of sym- 
metry is obviously some level crossing coming from the 
interplay of the current operator and the original Hamil- 
tonian. One might speculate that the occurance of such 
symmetry enhancements is not an accidental but a gen- 
eral feature of current carrying systems. 



B. Correlations 

The p^{r) correlations can be calculated easily and, 
just as in the equilibrium case, one finds /3^(r) ^ r~^. 
The difference from the equilibrium is that the oscilla- 
tory modulation of the decay (present in equilibrium 
for h ^ 0) becomes more complex. Such modulation has 
been observed in case of the imposed magnetic current 
(see Sec.|l|) as well as in the transverse Ising model with 
energy current ||^ . The exponent of the power law decay, 
however, is unchanged when the currents arc introduced 
in all of the above examples. Thus it seems that p^'{r) cor- 
relations are not too sensitive to the presence of currents. 
A possible reason for this apparent rigidity is, perhaps, 
the lack of internal interactions among the z-components 
of the spins. 

It is harder to calculate the p^{r) — p^{r) correlations 
but they show a more interesting behavior. Some of our 
results described below are exact and were derived by 
combining the Wick theorem and the spin rotation trans- 
formation (III) which relates the correlation functions be- 
tween states where the ground-state occupation pattern 
in fc-space is identical up to shifts fc — s- (k + a) mod 27r. 
As we shall see, these exact results are restricted to the 
boundaries of regions Q) and @ . At a general point 
{h,j^), we were able to calculate /9^(r) = p'-^{r;h,j^) 
numerically (for r < 100 lattice spacings) using the fact 
that the square of the correlation can be expressed as a 
determinant of a 2r x 2r matrix with exactly calculable 
elements. 

Let us start by enumerating the exact results. The 
boundaries of region (2) are discussed in points 1-3 while 
the boundaries of (T) are treated in points 3-5. 

1. As discussed in Sec.||, the correlations in the trans- 
verse XY model without current, i.e. on the line 
(0 < /i < 1, = 0), are known and the r — > oo 
asymptotics of p^ is given by 

p"(r; /i, j^-0)^C (1-/^2) 1/4 ^_i/2 _ (20) 

2. The correlations on the h = 1 line can be related to 
the equilibrium case (0 < h < = 0) and one 
finds: 



p^ (r; 1, j^) = r; V^l - n^j^\o] cos (|r) 



(21) 



Thus, the large-distance behavior is given by 



p^ (r; 1, f) = cV^r-'/^ cos(^r) . (22) 

Since the whole @ phase can be projected onto the 
h = 1 line, we find that correlations decay as r~^^'^ 
for h> 1. 
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3. The correlations on the 'high-symmetry' hne (i.e 
on the boundary between (I) and (g) ) can also be 
related to equilibrium: 

p''{r;h,f = h^/Ti) = p^(r;0,0)cos(fc,ir) 

- Cr-i/2cos(fchT) (23) 



so we find again a r ^1"^ decay in the r — > oo limit. 

4. On the line of maximal-current [ft, < i,j^iaj: ~ 
(1 -I- ft,^)/(27r)], the correlations can be expressed 
in terms of those on the line (h = 0, j^): 



P'^ir; h,j^^J =p'' [r; 0, I cos 



Unfortunately, this does not help in calculating the 
r — > oo behavior. 

5. The long range behavior of the correlations at the 
intersection point [0, l/(27r)] of the {h = 0,j^) and 
the (/i, Jmaj;) lines is also calculable: 







;0,o)]' 



J mteger 
otherwise 



(25) 



and it is remarkable that the correlation function 
in this point decays as 1/r instead of 

The exact results can be summarized as follows. The 
= p^ correlations decay as l/Vr on the boundaries 
of region @ while a 1/r decay can be observed in the 
upper-left corner [0, l/(27r)] of the phase diagram. 

Numerical calculations suggest, however, that the 1/r 
asymptotics is more general than it looks from the ex- 
act results. There is a strong indication that the large- 
distance asymptotics is actually 1/r everywhere in region 
(T) and @ apart from the boundaries of region @ . Fig.^ 
shows an example of numerical results at a general point 
of the phase diagram. The following formula gives an 
excellent fit to the numerical data throughout the phase 
diagram (except very close to the lines with the l/'s/r 
behavior): 



h=o.l j =0.38 

fit for r even 
fit for r odd 




FIG. 2. Absolute value of the p^{r) = (s|sj_^^) correlation 
function at a general point of the phase diagram. The distance 
r is measured in units of lattice spacing. The long-range de- 
cay of p^{r) is fitted by the expression ( p6| ) (solid lines) (for 
this point of the phase diagram one has ai ~ 02 ~ 0.0934, 
as ~ 0.0469). 



As one can see from (p6|), the amplitude of the 1/r 
decay is modulated with the critical wavenumbers kh 
and k\. On the 'high-symmetry' transition line we have 
kh = k\ and the transition across this line takes us from 
region ® where kh < k\ to region @ where kh > k\. 
Thus one can view the 'high-symmetry' line as a line of 
degeneracy where two characteristic wavelengths of the 
system become equal. 

This transition may resemble transitions arising from 
competing wavelengths but, actually, here we do not have 
a competition between kh and k\. Due to the product 
form of Afc , kh is independent of A while k\ is independent 
of h. Nevertheless, this transition does have similarities 
with second order transitions in that the correlations de- 
cay more slowly (1/r — > 1/Vr) and, furthermore, one can 
observe scaling upon approach of the transition line. In 
order to see this, let us assume that the distance from 
the transition line, kh — k\, provides the single diverging 
lengthscale which generates the 1/r ^ I/aA' crossover 
in correlations. Then one should observe scaling when 
plotting the following ratio: 



p^{r;h,j^) _^( \kh-kx\ ^ 



p-{r;K,j^) 



(27) 



xf\_f {o-i cos{khr) + a2 cos{kxr))/r for r even, 
^ 1 {a:jCOs{khr) + a2C0s{kxr))/r for r odd. 

(26) 

The Qi coefficients are functions of h and , and for the 
® phase ai = 02 seems to be valid. 



where (hcj^) is a point on the transition line and the 
{h hc,j^ Jc) limit is taken. Note that the long 
range behavior of the denominator {p^ on the phase tran- 
sition line) is known ( psf ) . As one can see from FigjI the 
data collapse is excellent thus supporting the assumption 
of scaling (|^. 
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scaling in the correlation function 



Jt/2 7C 

l(k,-k^)/2lr 



371/2 



FIG. 3. The absolute value of the scaled correlation 
function at four points of the phase diagram near the bound- 
ary between regions ® and @ . Two points are above and 
two are below this 'high'-symmetry line. The data points 
showing large deviations from scaling come from arguments 
where both the numerator and the denominator in (M) are 
close to zero. 

It is interesting to note that the scaling function ap- 
pears to be the same on the both sides of the phase- 
transition line. Furthermore, $ is independent of the 
crossing point {hcjj^) unless we are close to the zero- 
field equilibrium point {h — 0, — 0). In this sense we 
have the kind of universality usually observed in critical 
phase transitions. 

The equilibrium point {h = 0,j^ = 0) is the endpoint 
of the 'high'-symmetry line. At this point, the ground 
state symmetry is higher than on the line and so we may 
expect that, provided scaling was still present, the scal- 
ing function would be different. This is indeed what we 
observe. 



/9^(r; 0,0) 



(28) 



As shown in Fig.^, scaling is indeed seen and the scaling 
function is significantly different from that found away 
from the {h — 0,j^ — 0) point. 

The numerical results presented above (as well as other 
data gathered in our explorations of the phase diagram) 
suggest strongly that oc 1/r for generic current- 
carrying states. Slower decay, oc l/y/r is observed 
only on the boundaries of region (§) and the crossover be- 
tween the 1/r and l/\/r behaviors can be understood in 
terms of single-lengthscale scaling. It is intriguing that 
there is a simple correspondence between the types of 
decay of correlations and the "band structure" of the 
ground state. The lines of slower decay of correlations 
coincide with those lines where the ground state is build 
by a single band of excitations in momentum space, 
whereas in all regions of 1/r-decay, the filling pattern 
of the ground state splits into two separate bands (Fig. 
1). _ _ 

Regarding the interplay of currents and correlations 
these results leave us with the following conclusions. 
First, we find that the large-distance correlations are not 
necessarily increased by switching on a current. Second, 
it is found that the equilibrium power-law correlations are 
not destroyed by the current, only the exponent in the 
power law is increased. This strengthens previous ob- 
servations that currents and power-law correlations are 
intimately related. Third, we find that the increase of 
current may lead to interesting phase-transition like be- 
havior related to the increase of symmetry at special val- 
ues of the current. 
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FIG. 4. The absolute value of the scaled correlation 
functions on the h = line near = 0. 

Approaching the point {h = 0,j^ = 0) along the 
{h = 0,j^ 0) line, one has again a diverging length- 
scale proportional to l/k\ and one can search again for 



IV. FINAL REMARKS 

A general conclusion we can draw from the present 
study of the transverse XY model and from comparison 
with the results on the transverse Ising model |^ is that 
currents appear to generate and maintain power-law cor- 
relations. An interesting feature oi XY model which may 
also have some generality is the increase of the symme- 
try of the ground state at special values of the energy 
current. This feature should certainly be searched for 
in other models as well as in experiments. It should be 
recognized, however, that both the XY and the trans- 
verse Ising models are integrable and, consequently, they 
are special in that conductivity and, in particular, the 
thermal conductivity is ideal for them (not only at zero 
but also at nonzero temperatures) [[l5|-pT[ . Thus it is an 
important next step to find out whether nonintegrable 
models have the same connection between currents and 
power-law correlations and, furthermore, whether they 
show any additional general features. 
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